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Abstract: We study quantum electrodynamics with 2+1 dimensional massless Dirac 
fermion around a Coulomb impurity. Around a large charge with atomic number Z > 
137, the QED vacuum is expected to collapse due to the strong Coulombic force. While 
the relativistic quantum mechanics fails to make reliable predictions for the fate of the 
vacuum, the heavy ion collision experiment also does not give clear understanding of this 
system. Recently, the “atomic collapse” resonances were observed on graphene where an 
artificial nuclei can be made. In this paper, we present our nonperturbative study of the 
vacuum structure of the quasiparticles in graphene with a charge impurity which contains 
multi-body effect using bosonization method. 
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1 Introduction 


The quantum electrodynamics (QED) is the most precisely experimentally tested theory 
in today’s fundamental theories. In usual perturbation theory of QED, the expansion 
parameter is the fine structure constant a ~ However, in strong external field, since 
the interaction is correspondingly strong, the perturbation theory breaks down. From the 
result of relativistic quantum mechanics, the vacuum around an atom with large atomic 
number Z > 137 is expected to collapse. 

In non relativistic quantum mechanics, in the region where the potential is larger than 
the energy, the wave function falls off exponentially. Namely, all the incoming particles are 
reflected; that is, the reflection rate is R = 1 and transmision rate is T = 0. On the other 
hand, in relativistic quantum mechanics, when height of the potential Vo is larger than 
twice of particle mass 2 m, the reflection rate becomes larger than unity (R > 1), which is 
called the Klein tunneling[l, 2], This mechanism originates from the fact that the Dirac 
equation has both the positive and negative energy solutions as opposed to the Schrdingier 
equation. The same mechanism also prevents the electron to form bound states in a very 
strong attractive potential. In particular, an electron around a nuclei with a sufficiently 
large atomic number Z does not form a bound state due to the strong Coulomb potential, 
and fall into the nuclei. This phenomenon is called the atomic collapse, and has been known 
theoretically for a long time. However, since the atom with Z > 137 can be created for 
only a short time in heavy ion collision experiment, it is difficult to observe the phenomena 
experimentally at the quantitive level [3, 4], 

The situation has changed since the discovery of the graphene in 2004[5]. The electric 
structure of the graphene at low energy is known to be the same as that of the massless 
Dirac fermion. In addition, the effective fine structure constant a is about 300 times as 
large as that in the Quantum Electro Dynamics (QED). Due to this property, the essential 
point of the physics in the strongly coupled QED can be tested in the experiment using the 
graphene. Putting charged impurities on the graphene, one can realize a system similar to 
the large Z atom system, which enable us to observe the “atomic collapse”. 

The graphene is very thin, very light, very strong, and has very high electron conduc¬ 
tivity and made from carbon atoms which is a ubiquitous element on earth. Therefore the 
graphene is expected to serve as ideal device in future. In this point of view, understand¬ 
ing the response of the electron to the charged impurity in graphene is a very important 
problem and is studied actively. This system is well studied in one body quantum mechan¬ 
ics as the system of the two dimensional massless electron in Coulomb potential [6-8]. It 
is predicted that when the charge of the impurity exceeds a critical value Z cr , the wave 
function drastically changes. The massless fermion forms infinite number of quasibound 
states with negative energy, and the characteristic resonances appear in the local density 
of states (LDOS) of the electron[9]. Inspired by these theoretical studies, the scanning tun¬ 
neling microscope (STM) experiment was carried out and a characteristic peak in LDOS 
was measured[10]. 

However, the above theoretical studies do not take into account the many body effect 
which involve electron-positron pair creation. In the graphene case, since the pair creation 
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can occur with no cost of extra energy, the many body effect should not be neglected, which 
should be treated in the quantum held theory. Moreover, because of the large coupling 
perturbative approximation cannot be valid. Thus, this problem should be studied in some 
nonperturbative way. 

We analyze the held theory of 2+1 dimensional Dirac massless fermion around an 
external charge using the bosonization technique. In two dimensional theory, the fermion 
theory is converted to the boson theory[ll, 12]. It is known that a part of quantum effect 
of the fermion theory can be extracted from the classical boson theory. The bosonization 
method has been used to analyze the system with the fermion around monopole assuming 
that the classical boson theory captures the essential features of the quantum effect of the 
original fermion theory [13, 14]. The bosonization method is applied also to the atomic 
collapse problem in 3+1 dimensions [15]. We apply this method to the atomic collapse 
problem in 2+1 dimensions. 

Following the studies in 3+1 dimensions mentioned above, first restricting the gauge 
and the fermion held to s-wave held, we reduce the theory to 1+1 dimensional fermion 
effective theory. Next, we map the two dimensional fermion theory to the two dimensional 
boson theory. Then we solve the classical equation of motion for the boson held. As a 
result, we find the vacuum structure including the charge screening of the impurity charge. 

This paper is organized as follows. In section 2, the result of foregoing analysis in 
one body theory for the Coulomb impurity problem on graphene is briefly reviewed. In 
section 3, we will explain the s-wave approximation and the bosonization formalism pro¬ 
posed in Ref. [15]. In section 4, we will show the details about our study of vacuum 
solution and the result of our numerical analysis. Section 5 is devoted to summary and 
discussion. 

2 Review on Coulomb impurity on graphene 

In this section, we review the Coulomb impurity problem on graphene. The electronic 
properties of the graphene are described by the tight-binding model where interactions 
between different orbits are neglected. And it is assumed that the electron can hop to only 
the nearest neighbor site. In momentum space, the energy of electron becomes zero at two 
points (K and K'). The low energy effective theory is obtained by expanding the equation 
which the electron obeys around these points. It is known that the effective Hamiltonian 
takes the same form as that of massless Dirac fermion. That is, the fermionic low energy 
excitation obeys the Dirac equation 



( 2 . 1 ) 


and has the linear dispersion relation 



( 2 . 2 ) 


at low energy. The parameter vf in the above equation is the Fermi velocity which is 


roughly estimated as vf ~ m- Since vf plays the similar role as the speed of light c in 
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quantum electrodynamics, the effective fine structure constant for the fermionic excitations 
on graphene is a e ff ~ This means that the massless Dirac fermion on graphene is 

strongly coupled. 

The behavior of electron in a hydrogen like atom is studied in relativistic quantum 
mechanics. It is known that the bound state of electron and a point charge Ze cannot exist 
when Za > 1. For such a strongly coupled system, it is expected that the strong electric 
field makes the vacuum unstable since the strong Coulomb potential causes particle-hole 
pair creations. Such a phenomenon is called the “atomic collapse” and has been discussed 
for a long time. In the experimental side, the atomic collapse has been tested in heavy-ion 
collision. However the instability of large atomic number nuclei makes it difficult to observe 
the phenomenon clearly. 

In the graphene case, such a situation can be easily set up due to the large value of 
the effective coupling a ~ of the Dirac fermion. Recently, Wang and his collaborators 
studied the graphene system with Coulomb impurities with STM and observed the reso¬ 
nance like the quasibound state[10]. They put Ca dimers as impurity, and measured the 
local density of states (LDOS) of electron around the impurity. They showed that the peak 
appears in energy dependence of LDOS. The peak point is below the Dirac point when 5 
Ca dimers are put. According to them, this is the quasi-bound state expected in one body 
theory. The quasi-bound state spatially spread through about 10 nm around the center of 
Ca dimers in this experiment. 

In view of this STM experiment, it is now very important to study the graphene 
system with Coulomb impurities theoretically. In one body theory, the solution of the Dirac 
equation with Coulomb potential by a charged impurity can be exactly obtained[6, 7]. The 
behavior of the solution drastically changes when Za > 1/2. Because the electrons in 
graphene are massless, they do not seem to make bound state even in small Za. However, 
by introducing graphene lattice cutoff, the quasi-stable bound state is predicted to appear 
in strong coupling case. 

In Ref. [9], the existence of the quasi-bound state is semiclassically discussed. Here, 
we briefly review their discussion. The Hamiltonian for 2 dimensional massless fermion in 
Coulomb potential is 

H = a ■ p -—. (2.3) 

r 

When we write the square of momentum in terms of the radial momentum p r and the 
angular momentum j 

p 2 =p 2 r+j 2 /r 2 , (2.4) 


the Hamiltonian (2.3) leads to 


Pr = \ £ + 


Za + j 


£ + 


Za- j 


(2.5) 


where e is energy eigenvalue. The classically forbidden region where p% < 0 corresponds to 

_ Za- j Za + j _ 

r i = ——— < r < —j—j— = r 2 . (2.6) 
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= r 2 . 






Notice that if Za > j, there exist classically allowed region inside; that is, r < r\. Therefore 
in strongly coupled case, quasibound states can be found by imposing the Bohr-Sommerfeld 
quantization condition 



can be calculated. However, since the atomic collapse is a phenomenon which comes from 
pair creation effect, it should be analyzed in a way which contain nonperturbative multi 
body effects. In the following section, we will show the 2+1 dimensional massless fermion 
version of the bosonization formulation proposed in Ref. [15]. 

3 Approximation and Formalism 

In this section, we study the vacuum structure of the massless Dirac fermion system in 2+1 
dimensions around a Coulomb impurity. In order to analyze the system nonperturbatively, 
we employ the method proposed in Ref. [15] for the atomic collapse QED in 3+1 dimen¬ 
sions. We firs restrict the theory with s-wave electromagnetic field and the lowest partial 
wave electron field. Under this approximation, the theory is reduced to 1+1 dimensional 
effective theory with time and radial degrees of freedom. We then bosonize the effective 1 + 1 
dimensional fermion theory. Since it is known that the bosonized theory captures important 
part of the nonperturbative effect of the original fermion theory even at the classical level, 
we study the nonperturbative vacuum structure by constructing the classical solution of 
the bosonized theory. 

3.1 1+1D Effective Theory 

Since the gauge held is in 3+1 dimensions, we start from the following gauge action 



(3.1) 


where the charge density of impurity is spherically symmetric p(x ) = p(r, t ), and normalized 


as f d 3 xp(x ) = 1. The s-wave electromagnetic held takes following form 

A 0 (x) = a 0 (r,t), Ai(x) = ridi(r,t), 


(3.2) 


where f t = n/r is zth component of the unit vector in radial direction. In this approxima¬ 
tion, the gauge action becomes 



(3.3) 


When the graphene is on z = 0 surface and the electron is trapped on this surface, the 
action for the electron coupled with the gauge held is 



(3.4) 
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where ip is 2 component Weyl spinor. We take gamma matrices as 


7° = 0 - 3 , 7 1 = io 2 , 7 2 = -*cn- 


(3.5) 


Because we are considering z = 0 surface and using the s-wave approximation (3.2), y 3 
disappears from (3.4). From now on, i runs from 1 to 2. The fermion action becomes 


Sf = 


d 2 xdtip t [(i<9o + eao) + + ef^ai)] ip. 


(3.6) 


We expand the fermion field ip as 


Ip = -T= Y] v m ,a(r , #)T m *(<£>), 
V 7 " ^ 

v m,cr 

where cr = ±1 and m is half integer, and 

1 / e i(m-l/2)tp \ 

“ 1/tn \ ae^ m+1 / 2 ^ ) ' 

is normalized as 

d^p^f m id m m' • 

Using the relation 


O'if jUl m,a — &^m,c r, 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


we get 


@idi1p — <— ^ ( U (d r V m aij'i T ^m,- 

+ r V r 


(3.11) 


Therefore the action for fermion becomes 


m 


Sf= drdt ^2 [ v m,a {ido + ea 0 + a(id r + ear)} v m><T - iav^'^m-o 


(3.12) 


m,cr 


We restrict ourself to consider only the lowest (j = 1/2) partial wave, and define 1 + 1 
dimensional fermion 


'U j m. • — 


1 + i 1 — i 


+-U" 3 

2 2 


sign(m)u m _ 


Vm,+ j 

sign(m)iu m _ J ’ 


(3.13) 


where m = ±1/2. From now on, we take 


7° = cr 2, 7 1 = 7cr + 75 = 7°7 1 = o"3 
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(3.14) 






as 2 dimensional gamma matrices. Then we can rewrite 2 dimensional fermion action as 


S f = 


drdtt 

m=± 1/2 


{ 7 °(*^o + ea 0 ) + 7 1 («<9 r + ear)} 


2r 


(3.15) 


The last term represents centrifugal force. Unlike that of Ref. [15], we have a different 
coefficient of centrifugal force term and no mass term. 

We have to set the boundary condition for fermion field u m at r = 0 by requiring no 
singularity at r = 0. From (3.8), 


^ 1 / 2 ,+ — ^ 1 / 2 ,— ) ^- 1 / 2 ,+ + ^- 1 / 2 ,- (3.16) 

has ip dependence. If the coefficients of these are finite value at r = 0, the singularity arises. 
So, we set the boundary condition 


v m ,+( 0 , t) - sign (m)v m - ( 0 , t) = 0 . 


Written in 2D fermion u mi 


(1 - 7 °)tt m ( 0 , t) = 0 . 


(3.17) 


(3.18) 


On the other hand, since 


^ 1 / 2 ,++ ^ 1 / 2 ,-) ^- 1 / 2 ,+ — ^- 1 / 2 ,- (3.19) 

don’t have ip dependence, the coefficient of these can be finite at r = 0. Therefore we can 
also use the same boundary condition as Ref. [15]. 

By the way, in one body theory, the boundary condition is set not at r = 0, but at 
r = to[ 6 , 7], which is lattice cut off size of graphene. And the cut off plays very important 
role to discuss the drastic change of wave function and quasi-bound state in strong coupling 
region. In our case, however, even if we set the boundary condition at r = ro, we get the 
same result for r = 0. Therefore, here we set the boundary condition at r = 0. 


3.2 Bosonization 

We apply bosonization to this theory. Regarding interaction term as perturbation, we 
bosonize free fermion field to free boson field, 


u m (r, t ) 


/ _/M 1/2 f -i/V^exp [iy/Tr(<j> m (r,t) + <t> m (r,t))} \ 
V27t/ y N^expliy/ni-^m^t) + 4> m (r,t))\ )' 


(3.20) 


where 


cj)(x) = lim 

e ->0 


dse es <j)(s,t), 


(3.21) 


and N^ represents normal ordering at IR mass scale /r. From now on, we use the overdot and 
prime for time and spatial derivative, respectively. Because the action and the boundary 
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condition are almost the same as Ref. [15], we can bosonize this theory following the same 
calculation. 

In this case, we should impose the boundary condition on the boson field. The boundary 
condition (3.18) is rewritten in boson field as 

</>m(0,t) = 0. (3.22) 


Free boson held can be expanded in plane wave as 


0(®, t) = f — \d\k)e ik ^ x+t) + a(fe)e ifc M + a{k)e ~ ik < x+t ) + a\k)e~ ik ^ 

Jk>o 27T 1 


(3.23) 


where a, a, at, a t are creation-annihilation operators satisfying appropriate commutation 
relations. While a(k), a(k) are independent operators without the boundary condition, with 
the boundary condition (3.22) 


0 = ^( 0 ,*) 


dk 


Jk >0 2 ^" 

these are dependent on each other 


(a^(fc) + a) (k))e lkt + ( a(k ) + a(k))e 


—ikt 


(3.24) 


a{k ) = — a(k). 

Then the boson held <f> and 0 can be written as 

dk 


0(r, t) = / — a(k)(e ikr - e- ikr )e~ ikt + a\k)(e~ ikr - e ikr Y kt 

Jk >o 2?r L 


and 


0(r, t) = [ —\(e ikr + e~ ikr ) a(k)e~ ikt + ( e ikr + e~ ikr ) a\k)e ikt 

Jk >o 2?r LV ' 


We split these into 


0 ( +)(r,t)= [ —a(k)(e ikr -e- ikr )e~ ikt , 
Jk> o 2vr \ J 

<t>(~\r, t )= [ —a\k)(e~ ikr -e ikr )e ik \ 
Jk> o 2vr V J 

0(+)(r,*) = [ —a(k ) (e ifcr + e" iH , 

A->o 2tt V / 

0(-)(r,t)= [ —a\k)(e ikr + e~ ikr )e ikt . 
Jk> o 2vr V / 


(3.25) 

(3.26) 

(3.27) 

(3.28) 

(3.29) 

(3.30) 
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(3.31) 










From the commutation relation [a(k),a^(k)] = 2n-^5{k — k'), we get the relation 

[4>^ (r, t) + 7(r, t), cf>^ (/, f') + rj'4 (/, f')] 

= - ^X 1 - X)^+ + (1 + h)(l + v')A- + (1 - ??)(1 + r/)B+ + (1 + rj)( 1 - r/)5_], 


where 


A±(r, t; r', t') = — [ dk^e ik ^ r - r '^ t ~ t '^ 


' fc >0 


= lim In (i/i[t — t! ± (r — r 7 ) — ie \), 


(3.32) 


(3.33) 


B±(r, t; r', t') = - I dk\e 


/ fc >0 


= lim In (i/i[f — ± (r + r') — ze]), 


(3.34) 


are renormalized at IR mass scale /i. B + _ arise from the boundary condition. 

Using the commutation relation (3.32), we rewrite the interaction terms in fermion 
theory in terms of boson held. After some point splitting procedure, we get 


u^u = — ^=e^d u (j), 
Wn 


(3.35) 


177577 = --—Nn cos(2v / vr</>), 

Zirr 


(3.36) 


where e is anti-symmetric symbol with e 10 = 1. 

We rewrite the fermion action (3.15) in terms of the above boson operators. In ai = 0 
gauge, 


/ ” 1 G- 1 

drdt ^2 -^=OQ0 m + cos(2 v / 7r0 m ) 

m=± 1/2 L y/TT nr 

where the second term is integrated by parts. Therefore we get the total action 


(3.37) 


S= drdt 


2tt rV 0 2 + e<F(r, t)a' 0 + ^ ~ A^a'o4>m + ^2 cos(2 y/n(f>m)j 

(3.38) 

where the 4>(r, t) is defined by 

<b'(r, t) = 4nZr 2 p(r, t). (3.39) 

From the action (3.38), we notice that ao has no dynamical degrees of freedom. Using 
the equation of motion for Oq 

4nr 2 a' 0 + e$(r, t) - -^=<j> m = 0, 

^' \ n 

m=± 1/2 v 
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we can eliminate ao- Therefore the Hamiltonian becomes 



we shift the energy so that the energy becomes zero when i j> m = 0 which is vacuum con¬ 
figuration with Z = 0. In the next section, we numerically calculate the solution which 
minimize this Hamiltonian. 


4 Study of Vacuum Solution 


In this section we find the classical solution which minimizes the bosonized Hamiltonian 
in the previous section. For this purpose, we have to solve the Euler-Lagrange equations 
for boson fields under the appropriate boundary conditions. The boundary condition at 
r = 0 is determined by eq.(3.22). The solution is characterized by the boundary condition 
at r = oo. 

Eq.(3.35) indicates that the density of electron p e (r) can be written in terms of the 
boson field as 

Pe(r) = = V] -)=d r ct>m{r). (4.1) 

Therefore, we get the spatial distribution of induced electron density corresponding to the 
solution. Total induced charge which screens the impurity charge is given by 


Qem = —e 


drp e (r ) 



5>m(oo). 

m 


(4.2) 


In order to study the vacuum structure, we consider only static solution 7r m = 0. We rewrite 
the Hamiltonian in terms of the new variable 


<t>± 


± 0 — 1 / 2 )- 


(4.3) 


Using the formula 

cos [v / 2vr(0+ + 4>-)} + cos [\/27r((/>_)_ - 4>-)\ = 2 cos (y / 27T(/>+) cos (\/2v T(j)-), (4.4) 










the Hamiltonian becomes 


where a 


^0+ + 4 >~) + 7^2 U _ cos(\/2vn/>+) cos(V2vr^-)} 

2 

|^. The Euler-Lagrange equations for </> + ,0_ are given by 

0 + - sin(v / 27 r</>+) cos(\/ 2 vr^_) - ^ = 0 , 



(4.5) 


(4.6) 


4>'!_ - - cos(^+) sin(v / 27T0_) = 0, (4.7) 

\/2iTr 2 

respectively. Since it satisfies Eq.(4.7) we can take the symmetric ansatz = 0. Then 
Eq.(4.6) reduces to 


K ~ sh sin(v/ ^ +) - d (*+ ~ = a 

We assume that the impurity charge is spherically spread over radius R: 

The corresponding <L(r) is 

H r) = { Z ^f ^<R) 

W \ Z (r> R). 


(4.8) 


(4.9) 


(4.10) 


Since Eq.(4.8) is a second order differential equation, in addition to the boundary condition 
at the origin we need to impose another boundary condition at r = oo. For finiteness 
of total energy, the boson held should asymptotically be constant (</>+ —> 0*) at large r. 
Substituting = q i>* into the Euler-Lagrange equation at large r, we find that 0* should 
be the solution of the following equation: 


sin 


(V2tt4>*) = —2a 



(4.11) 


Notice that the asymptotic value y f 0* can take non-integer value. Charge screening with 
non-integer charge may seem counter intuitive if one tries to interpret the phenomena as 
particle hole pair creation. One should interpret such screening as the polarization effect. 

In fact, it is known that the screening of non-integer charge actually occurs in massless 
Schwinger model [13, 16, 17]. In the following subsections, we show the detailed numerical 
analysis and its results. 
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4.1 Numerical Analysis 

4.1.1 Strategy 

Our numerical analysis is done in various parameters a, Z. according to the following steps: 

(i) Find the solution of Eq.(4.11) and obtain the asymptotic form at large r. 

(ii) Solve the Euler-Lagrange equation (4.8) with the boundary condition at large r (4.17) 

with various A. 

(iii) Find A with which the solution satisfies the boundary condition at r = 0 (3.22). 


4.1.2 Asymptotic form 

In order to numerically solve the Euler-Lagrange equation (4.8), we should find the asymp¬ 
totic form at large r. To do so, we parameterize (j>+(r) by introducing a function / which 
describes the deviation of (j>+(r) and 0* at large r as 

<M r ) = <t>* - fir). (4.12) 


where 0* is the solution of Eq.(4.11). Substituting Eq.(4.12) into Eq.(4.8), and expanding 
it up to linear order in /, we obtain 

f" ~ ~2 (^os(V2ncj)*) + / + 0(/ 2 ) = 0 (4.13) 

Assuming that the solution for / can take the form 

/( r) « (4.14) 

at large r with A being a constant and substituting it into Eq.(4.13), we find that the power 
A satisfies the following equation: 


A 2 + A- 



= 0 . 


From Eq.(4.15), A should be 


A 


1 

2 



(4.15) 


(4.16) 


We show some examples of Z = 4 case. In this case, there are three screening patterns 
depending on the value of a as shown in Fig. 1. In a < 0.14 case, there are five values of 
</>*. However, when 0* is equal to the second or fourth smallest value, based on Eq.(4.16), 
A becomes imaginary. Only the solutions with the real positive values of A make sense. So, 
there are three possibilities. For 0.14 < a < 0.34 case, there are three values of 0*. Similarly 
the second smallest value of is not a physical solution. So, there are two possibilities. 
And in 0.34 < a case, there is only value for 0* which corresponds to the full screening 
solution. 

In Fig. 2, we show the number of possible asymptotic solutions at large r for each set 
of values of (a, Z). 
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Figure 1. The green line is l.h.s. of eq.(4.11), and blue, red, yellow lines are r.h.s. of eq.(4.11) 
with a = 0.1, a = 0.2, a = 0.4, respectively. 

4.1.3 Example of the solution 

Starting from the asymptotic solutions and solving the differential equation numerically, 
we can obtain the full solution. Taking the following asymptotic from 

4>+{r) « </>* - (4.17) 

at large r and varying A, we can search for the physical solution which satisfies the boundary 
condition at r = 0. Practically, we solve Eq. (4.8) from r = 0.0017? to r = 1000007?, 
setting the boundary condition at large r = 1000007? with various values of a, Z. For 
illustration, we show the example for Z = 4 and a = 0.2. In this case, there are two 
asymptotic solutions, but only the solution which realizes the smallest value of </>* can 
satisfy appropriate boundary condition. The full solutions from the other asymptotic forms 
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Figure 2. Number of possible asymptotic solutions at large r for each set of values of (a, Z). 
Crosses are the parameter points where we solved Eq. (4.8) 
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do not satisfy the boundary condition at r = 0 but end up have positive values no matter 
how we choose the value of A. We show the solution of < j>+(r) in Fig. 3. In the other case, 
the shapes of solutions are qualitatively similar to the solution in this case. The induced 
electron density is depicted in Fig. 4. We notice that most of the induced electrons fall 
into the impurity. 



Figure 3. The solution for a = 0.2, Z = 4. 



Figure 4. Induced electron density with a = 0.1 
(blue line), a = 0.2 (red line), a = 0.4 (yellow 
line). 


4.2 Result 

4.2.1 Phase structure 

We looked for the solution for various set of parameters of (a, Z ), where the parameter set 
is given in Fig. 2 . We found that only the solution with the smallest value of can satisfy 
correct boundary condition at r = 0 in all cases. From this fact, we reach the conjecture that 
the magnitude of screening can be determined by the smallest intersection of sin(v / 27^|?^ , *) 
and — 2a{\j2/'K(f) Jf — Z). According to this conjecture, we get effective impurity charge seen 
from infinitely separated point, 


Z e s = Z - 



(4.18) 


which is screened by induced charge Fig. 5. Notice that when a > 0.2, the effective impurity 
charge Z e g in any odd Z case is the same one as in Z = 1 case. Also when a > 0.14, Z e g 
in any even Z case is the same one as in Z = 2 case. From this result, a phase diagram 
of screening is described as in Fig. 6. In larger Z case, more branches appear in small a 
regime. 


4.2.2 Scaling law 

The induced 2D electron density can be obtained as 


n[r) 


Pe(r) 

2irr 


(4.19) 
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Figure 5. Z e s for each Z. Dashed lines in Z = 3,4 cases describe Z e g in Z = 1,2 cases, respectively. 



Figure 6. Phase diagram of screening. Crosses are the parameter points where we solved Eq. (4.8) 


If graphene sheet can be treated as perfect metal, the scaling law is calculated as in Ref. 

[ 18 ]: 


n(r) oc r 


(4.20) 
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in the range of distances 1 <C r/R <C 2 a 2 Z. In our calculation, we fit the scaling law 


n(r) oc r 7 


(4.21) 


in the range of distances 1 <C r/R <C 10. The scaling exponent 7 depends on parameters 
a, Z as shown in Fig. 7. In small screening regime, we get 7 ~ 2.7, independently of 
a. Near the value of a where magnitude of screening jumps, 7 drastically decreases. In 
larger a regime, 7 increases and becomes close to the value calculated in perfect metal 
approximation. 


7 


7 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 


a 


a 


Figure 7. Crosses describe the scaling exponent in each Z case. Dashed line describes one in 
perfect metal approximation. 


5 Summary and Discussion 

In this paper, we studied quantum held theory with the 2+1 dimensional massless fermion 
around an external Coulomb held. We reduced the theory to a two dimensional fermion 
theory, where the higher partial waves are neglected. Bosonizing the theory, we have found 
the static solution of classical equation of motion for the boson held. The magnitude 
of screening is determined only by the asymptotic equation of motion. Which of these 
asymptotic solutions satishes the boundary condition at r = 0 is determined by dynamics. 

Through the study of several examples, we have concluded that the realized solution 
is always the smallest screening one. As a result, we have found patterns of screening 
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depending on the coupling a and the impurity charge Z. The screening charge undergoes 
a drastic change as we change the value of a at some critical values. We also obtained the 
phase diagram characterized by the patterns of screening. 

By solving the equation of motion in full spatial regime, we have obtained the spatial 
distribution of density of the induced electron. The radial profile of the two dimensional 
induced charge density can be fitted by negative power in r which is the distance from the 
impurity. In weak coupling regime, scaling exponent 7 is independent of a and Z\ 7 ~ 2.7. 
Near the screening jumping point, 7 decreases. This means that the induced fermion is 
widely spread near the screening jumping point. And in larger a regime, 7 become close to 
the value of the perfect metal approximation; 7 ~ 3. 

The validity of the approximation to neglect higher partial wave can be discussed 
somewhat in semi classical theory mentioned in section 2. According to the semi classical 
theory, only Za > j wave can form quasi-bound states. So, the fermion mode whose 
angular momentum j is higher than Za is irrelevant to anomalous behavior of the electron 
in strong Coulomb potential. When Za > 3/2, the next to lowest partial wave j = 3/2 
should be relevant to this problem. Therefore our approximation should be valid only when 
Za < 3/2. 

To compare our analysis with the result of one particle theory or the experiment, 
many things remain to be done. Validity of classical treatment for boson theory should 
be confirmed quantitatively. In Ref. [19], the bosonized atomic collapse problem in 3+1 
dimensions is treated within small fluctuation approximation. They show the existence of 
meta stable states in supercritical phase. In the same way it may be possible to show the 
existence of the meta stable states in our 2+1 dimensional massless fermion case. 

The contribution of higher momentum partial wave should be evaluated for under¬ 
standing larger a, Z case. Furthermore to understand the behavior in the regime closer to 
the impurity, the effect of graphene lattice should be considered. For that purpose, the 
simulation by lattice gauge theory is important. 
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